Low density approach to the Kondo-lattice model 
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We propose a new approach to the (ferromagnetic) Kondo-lattice model in the low density region, 
where the model is thought to give a reasonable frame work for manganites with perovskite struc- 
ture exhibiting the "colossal magnetoresistance" -effect. Results for the temperature- dependent 
quasiparticle density of states are presented. Typical features can be interpreted in terms of elemen- 
tary spin-exchange processes between itinerant conduction electrons and localized moments. The 
approach is exact in the zero bandwidth limit for all temperatures and at T = for arbitrary band- 
widths, fulfills exact high-energy expansions and reproduces correctly second order perturbation 
theory in the exchange coupling. 
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I. INTRODUCTION 

The Kondo-lattice model (KLM) [1] describes the in- 
terplay of itinerant electrons in a partially filled en- 
ergy band with quantum mechanical spins (magnetic mo- 
ments), localized at certain lattice sites. Characteristic 
model properties result from an interband exchange in- 
teraction between the two subsystems. 

On the one hand, the energy bandstructure is modi- 
fied by the magnetic state of the spin system (temper- 
ature dependences, band splittings, band deformations), 
while, on the other, the magnetic state of the spin system 
is even provoked by the itinerant electrons because the 
KLM does not incorporate a direct exchange between the 
moments. The model-Hamiltonian consists of two parts 



Hs is the kinetic energy of itinerant band electrons. 



H, 



Y.T. 



(1) 



(2) 



I'ja 



where c\^{ci„) is the creation (annihilation) operator of a 
band electron specified by the lower indices. Tij are the 
hopping integrals. The second term in (1) is an interband 
exchange term with coupling strength J, written as an 
intra-atomic interaction between the conduction electron 
spin Gi and the localized magnetic moment represented 
by the spin operator S^: 



Hsf = ^J/_^'^i ■ Sj 



(3) 



According to the sign of the exchange coupling J, a par- 
allel (J > 0) or an antiparallel (J < 0) alignment of itin- 
erant and localized spin is favoured with remarkable dif- 
ferences in the physical properties. The parallel (J > 0) 
orientation is often referred to as " ferromagnetic Kondo- 
lattice model" (FKLM); alternatively known as "s-f or 
s-d model". The applications of the KLM are rather 
manifold. 



A. Magnetic semiconductors: 

Prototypes are the Europium chalcogenides EuX (X 
= O, S, Se, Te) [2] which are known to exhibit a spec- 
tacular temperature dependence of the band states. The 
"red shift" of the optical absorption edge upon coohng 
from T — Tc to T = OK [2, 3] is due to a corresponding 
shift of the lower conduction band edge. There is clear 
evidence that in these materials the exchange J is pos- 
itive, typically of order some tenth of eV. The coupling 
can therefore be classified as weak to intermediate. 



B. Semimagnetic semiconductors: 

In systems like Cdi-xMuxTe and Hgi-^FexSe ran- 
domly distributed Mn'^'^ or Fe^~^ ions provide local- 
ized magnetic moments which influence, via the exchange 
mechanism J, the band states of the II- VI semiconductors 
CdTe and HgSe. For moderate doping x, the moments do 
not order collectively so that a striking temperature de- 
pendence, as that of the magnetic semiconductors (EuX), 
cannot be expected. However, an anomalous magnetic 
field dependence of optical transitions and therewith of 
the bandstructure is observed [4] ("giant Zeeman split- 
ting"). From respective experimental data, J > can be 
concluded. The coupling must be classified as weak. 



C. Local-moment metals: 

In ferromagnetic metals such as the Rare Earth el- 
ement Gd, the magnetism is due to strictly localized 
4f electrons while the conductivity properties are deter- 
mined by itinerant (5d, 6s) electrons. The T — 0-moment 
of Gd is found to be 7.63/iB [5]. 7/is stem from the ex- 
actly half-filled 4f shell. The excess-moment of 0.63/iB 
originates from an induced spin polarization of the " a pri- 
ori" non-magnetic conduction bands, indicating a weak 
or intermediate coupling J > [6] . Many of the recent 



research activities have been focussed on the tempera- 
ture dependence of the induced exchange sphtting. Is 
it collapsing for T ^ Tc or does it persist even in the 
paramagnetic phase [6, 7]? The J- induced correlation 
and quasi particle effects in the valence and conduction 
bands of Gd (or equivalently Dy, Tb) lead to highly com- 
plex and therefore controversially discussed photoemis- 
sion data [8, 9], the interpretation of which is far from 
settled (see review in [7]). While the magnetic ordering 
of the semiconductors and insulators (class 1) has to be 
explained via special supcrcxchange mechanisms, which 
is beyond the field of application of the FKLM, it is com- 
monly accepted that the collective magnetism of the " lo- 
cal moment" metals is caused by the RKKY interaction. 
The latter is also based on the exchange interaction J. 
The FKLM therefore provides, at least in a qualitative 
manner, a self consistent description of magnetic and elec- 
tronic properties of materials such as Gd [6, 10]. 



D. Manganites-perovskites 

Since the discovery of the " colossal magnetoresistance 
(CMR)" [11, 12], the manganese oxides with perovskite 
structures Ti^^D^MnOz (T = La, Pr, Nd ; D = Sr, 
Ca, Ba, Pb) have attracted high scientific interest. The 
prototypes Lai-x{Ca, Sr)xMn03 are since long the pro- 
tagonists of the "double exchange" mechanism [13]. Re- 
placing, in La^^Mn^^Os, a trivalent La"^+ ion by a diva- 
lent earth-alkali ion (Ca^+,S'r^+) requires an additional 
electron from the manganese for the binding. The result 
is a homogenuous valence mixture of the manganese ion 
{MnltxMni+) . The three 5d - t2g electrons of Mn'^+ 
are considered as more or less localized forming a local 
S — 3/2 spin. The fourth electron in Mn^^ is of 'Sd — Cg 
type and is itinerant. It is assumed that it interacts 
via intra-shcU Hund's rule coupling ("double exchange 
model" [14]) with the S ~ 3/2 spins. The manganites are 
bad electrical conductors. It has therefore to be assumed 
that the intraatomic coupling J > is much stronger 
than the hopping-matrix element |t| (J >> |t|). Theo- 
retical estimates for the bandwidth yield W = 1 — 2 eV 
[15-17], experimental data propose W = 3 — A eV [18, 
19]. The exchange coupling J is not very well known, 
the J = 1 eV oi refs. [15, 20] are sometimes questioned 
as being too small [21]. In any case, the manganites 
belong to the strongly coupled FKLM which cannot be 
treated perturbatively with respect to J. The FKLM will 
certainly be overcharged to reproduce all the details of 
the rich phase diagram of Lai-xCaxMnO^, e.g., accord- 
ing to which the ground state is antiferromagnetic for 
a; = and 1, ferromagnetic for x « 0.2 — 0.4, with para- 
magnetic regions and phase separations in between [12]. 
Nevertheless, the FKLM is thought to give a reasonable 
frame work for an at least qualitative understanding of 
the interesting physics of the manganites [22, 23] 



E. Heavy Fermions 

The above subclasses are all characterized by a fer- 
romagnetic exchange interaction J > 0. The original 
Kondo-lattice model [24], however, refers to J < 0, 
favouring an antiparallel alignment of conduction elec- 
tron spin and localized spin. This situation is obviously 
realized in the Heavy-Fermion systems, which are to be 
found especially among Ce-compounds and which have 
provoked intensive research activities because of their ex- 
traordinary physical properties. Doniach [24] was the 
first to point out that there should be a phase transition 
from a magnetic state for small \J\ to a non-magnetic 
Kondo state for large \J\ characterized by a screening of 
the local moments by the conduction electron spins. The 
magnetic state is due to the RKKY-interaction, which as 
an effect of second order (~ J^) is independent of the sign 
of J. However, the Kondo screening is of course absent 
for J > 0, i.e., for all the above discussed subclasses. For 
most of the Heavy-Fermion systems, the RKKY coupling 
favours an antiferromagnetic ordering of the local mo- 
ments. In CeCuQ^xAux the competitive behaviour of the 
RKKY and the Kondo screening tendencies can impres- 
sively be observed by varying the concentration x [25]. 
CeCue{x = 0) is non-magnetic because of perfect Kondo 
screening, while for a; > 0.1 the RKKY component dom- 
inates taking care for an antiferromagnetic ordering up 
to a; = 1 {CeCu^Au) with increasing Neel-tempcrature 
Tn for increasing x. 

J < does not necessarily lead to antiferromagnetism. 
The compound CeSix is ferromagnetic for 1.6 < a; < 
1.85 [26] with a strongly reduced magnetic moment. The 
Curie temperature of the ferromagnetic Kondo system 
(J < 0) CeNixPti-x first increases in between a; = 
and X = 0.5 from 5.8K (a; = 0) to about 8.6K (a; = 
0.5) in order to decrease then rapidly and disappearing 
eventually at a; = 0.55 [27]. The magnetic moment per 
Ce ion diminishes steadily with x because of increasing 
Kondo screening and disappears completely at a; = 0.95. 

The above-presented list documents the rich variety of 
applications for the KLM. Since the many-body problem 
of the Hamiltonian (1) could not be solved exactly up 
to now, approximations must be tolerated. Most of the 
recent theoretical papers, aiming at the CMR-materials, 
assume classical spins S ^ oo [28-30] , mainly in order to 
be able to apply " dynamical mean field theory" (DMFT) 
to the FKLM problem. The merits of DMFT, e. g., with 
respect to the Hubbard model, are indisputable, but the 
assumption of classical spins in the KLM appears very 
problematic. Several important features, as, e.g., the 
magnon emission and absorption by the itinerant elec- 
trons, are excluded from the very beginning. The impor- 
tance of such effects has been discussed in detail in rcf. 
[10]. Conclusions such as, that at T = the spins of the 
Cg electrons are oriented parallel to the t2g— spins [30], 
are correct only for 5* ^ oo. For any finite spin, there 
is a considerable amount of | —spectral weight overlap- 
ping with I —states even for very large J. Recently, a 



DMFT-based approach to the KLM with quantum spins 
has been proposed [31] which uses a ferniionization of the 
local spin operators. The theory is restricted to 5* = 1/2 
but retains the quantum nature of the spins. A band 
splitting, which occurs already for relatively low inter- 
action strengths, can be related to distinct elementary 
excitations, namely magnon emission and absorption by 
the itinerant electron and the formation of magnetic po- 
larons. The results, which are in remarkable agreement 
with those from the "moment conserving decoupling ap- 
proach" (MCDA) in ref. [10], confirm the importance of 
the quantum nature of the spins. 

Due to some reasons, the above-mentioned theories 
[10, 31] are best justified for weak and intermediate cou- 
plings J. In this paper, we propose an approximate 
scheme which mainly aims at the strong coupling regime 
{JS >> W : W is the bandwidth) being nevertheless 
perturbationally correct up to order J^. The idea is to 
construct a sclfencrgy ansatz which interpolates between 
exactly known limiting cases and reproduces the correct 
high-energy expansion of the selfenergy. To demonstrate 
the method as clearly as possible, we restrict our consid- 
erations to the low-concentration region, performing the 
detailed calculation for a single electron in an otherwise 
empty conduction band. The theory is outlined in Sect. 
II, while Sect. Ill brings a discussion of the results. 



II. THEORY 

A. The many body problem 

The model-Hamiltonian (1) defines a non-trivial many- 
body problem, the exact solution of which is known only 
for a small number of special cases. For practical rea- 
sons, it is sometimes more convenient to use the second 
quantized form of the exchange interaction (3): 
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(4) 



is of primary interest. Its equation of motion reads: 

m 

(7) 

where the two types of interaction terms in (4) lead to 
the " spinflip function " , 



F,„,^,,iE) = ((Src. 



JO- 



and the " Ising function" : 



(8) 



(9) 



The two "higher" Green functions on the right-hand 
side of (7) prevent a direct solution of the equation of 
motion. A formal solution for the Fourier-transformed 
single-electron Green function. 



G'k.(^) = ((ck.;c. 



kcr 



E S-e(k)-Ska(S) 



(10) 



defines the in general complex self-energy Sko- (E) by the 
ansatz 

Hsf, ck.]_ ; 4,))^ - Sk,(S)Gk.(^). (11) 

e(k) are the Bloch energies: 

An illustrative quantity which we are going to discuss 
in the following is the quasiparticle density of states (Q- 
DOS): 



'Ae)^--^J2^^g^Ae+'^' 



(13) 



Here we have used the abbreviations: 



For the general case neither Eko-(-£') nor G]ia{E) can be 
determined exactly. However, some rigorous statements 
are possible and shall now be listed up. 
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0(j^ 5 ^j — '-'7 ' 'i^a'^j • y^) 



B. Zero- bandwidth Hmit 



The first term in (4) describes an Ising-like interaction 
between the z-components of the localized and the itin- 
erant spins. The second term refers to spin exchange 
processes between the two subsystems. 

If we are mainly interested in the conduction electron 
properties, then the single-electron Green function. 



Gija{E) - 



(6) 



The final goal of our study is to arrive at a self- 
energy formula being credible first of all in the strong 
coupling limit {JS » W). That means, in particular, 
that our approach has to fulfill the exactly solvable zero- 
bandwidth case [32]: 



Tij -^ TaSij ; e(k) -^ To Vk. 



(14) 



The conduction band is shrunk to an N-fold degenerate 
level Tq. The localized spin system, however, is furtheron 



considered as collectively ordered for T < Tc by any di- 
rect or indirect exchange interaction. The latter is not 
a part of the KLM. The localized magnetization (S^) 
Therefore enters the calculation as external parameter. 
With (14), the hierarchy of equations of motion for the 
single-electron Green function Gij (E) , following from eq. 
(7), decouples exactly [32]. The result is a four-pole func- 
tion: 
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with spin-independent poles at 

Eia =To- ^JS ; E2a =To + i J(^ + 1). 



(15) 



(16) 



^3. =To- i J(5 + 1) ; Ei, = To + ^JS (17) 

The j = 1, 2-excitations (Eq. (16)) refer to singly occu- 
pied sites; more strictly, they appear when the test elec- 
tron is brought to a site, where no other "conduction" 
electron is present. It then orients its spin parallel (Ei^) 
or antiparallel (i?2cr) to the local spin. These excitations 
are bound to spin-dependent spectral weights: 

-(5+l)(n_,)} 



1 



a2cy = 



25 + 1 
Here we have abbreviated: 



{5-m, -A_, -5(n_,).} (19) 



(20) 



The expectation values in the spectral weights ai^2a 
are, therefore, all selfconsistently determinable by the re- 
quired single-electron Green function itself. 

The two other poles E^^ andi<^4CT are bound to double 
occupancies of the lattice site. The test electron enters a 
site which is already occupied by another electron with 
opposite spin. The corresponding spectral weights. 
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vanish in the limit of zero band occupation. It may be 
considered a shortcoming of the KLM that the excitation 
energies (17) do not contain the Coulomb interaction en- 
ergy. Switching on a Hubbard interaction U leads to an 
additive term U in E^^ as well as in E/^a- [32], shifting 
these excitations to higher energies. While the Hubbard- 
U is of course the exact ansatz in the zero-bandwidth 
limit, it is not so obvious by which type of Coulomb inter- 
action the KLM should be extended ("correlated KLM" 
[30]) when aiming at one of the subclasses described in 
the Introduction. To avoid this ambiguity we restrict our 
following considerations to the low-density limit {n -^ 0), 
where the selfenergy of the zero-bandwidth KLM reads 
according to eqs. (16)-(19): 
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J^SiS+l)~^Jm,{E-To) 
E-n- \j{m, + 1) 



(26) 



This rigorous result will be exploited later for testing our 
approximate theory. 



C. Ferromagnetically saturated semiconductor 



y^i ^i-a^i<y) + ^o- \pi n-ia) 



(21) 



The " mixed" correlation function A^- can be derived via 
the spectral theorem from the Ising- and the spinflip- 
functions (8 and 9). Exploiting the equation of motion 
(7), this can even be expressed in terms of the single- 
electron Green function: 

1 1 f^°° 

A. = -^^5]y dE,UE)iE-e{k))ImG^,{E) 

(22) 

where f-{E) = (1 + e'^"''))"^ is the Fermi function (^ 
is the chemical potential). Similarly it holds for the spin- 
dependent particle numbers: 

1 1 f^°^ 

M = -^j^Y.j_ dEf^{E)ImG^AE). (23) 



There is another very instructive limiting case that can 
be treated exactly. It concerns a single electron in an oth- 
erwise empty conduction band interacting with a ferro- 
magnetically saturated local moment system (T = 0). In 
the zero-bandwidth limit (Sect. II. 2) for the j-spectrum, 
all the spectral weights (19), (24) and (25) disappear, 
except for aij = 1. In the ^.-spectrum the levels Ei^ 
and E21 survive with the weights aii — 25+T ^^'^ 

•^24 — 25+1 ■ 

For finite bandwidth, the mentioned special case is that 
of a ferromagnetically saturated semiconductor (EuO at 
T = 0!) [10, 31, 33-35]. In this situation, an t-clcctron 
has no chance for a spin-flip, the corresponding quasi- 
particle density of states, p-f (_E), is therefore only rigidly 
shifted compared to the "free" DOS [10] and the self- 
energy is a constant: 



^(T=0,n=0) 
■"kT 



_ „(T=0,n=0) 



(E) = ^- 



{E) = -ljS. 



(27) 



The j-spectrum is more complicated since a j-electron 
has several possibilities to exchange its spin with the an- 
tiparallel, localized spins. The spinflip function (8) does 
not at all vanish as in the f-case. Nevertheless, the prob- 
lem is exactly solvable resulting in a wave-vector inde- 
pendent self-energy: 



^{T=0,n=a)i 
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Go{E) is the "free" propagator: 



Go 
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(28) 



(29) 



The reason for the wave-vector independence of the 
self-energy can be traced back [10] to the lack of a direct 
(Heisenberg) exchange term in the model-Hamiltonian 

(1). Therefore Y,\ ~ '"~ {E) does not contain magnon 
energies ?ia;(q) which come into play when the excited 
J,-electron flips its spin by magnon emission. Neglecting 
the exchange between the local-moment spins Si may be 
considered as the ^^^(q) = case. As a consequence, the 
electronic self-energy becomes wave-vector independent. 
There does not arise any problem in calculating the limit 
(n = 0, T = 0) with the inclusion of a Heisenberg ex- 
change (~ JijSi ■ Sj). Then the wave- vector dependence 
of the selfenergy reappears [10]. 



D. Second-order perturbation theory 

Conventional diagrammatic perturbation theory for 
the Kondo-lattice model does not work because of the 
lack of Wick's theorem. A fertile alternative is the Mori- 
formalism [36, 37], which allows for a systematic expan- 
sions of the electronic selfenergy of the KLM with respect 
to the powers of J. That has successfully been done pre- 
viously for the weakly coupled Hubbard model by the 
use of the modified perturbation theory of [38, 39]. In 
the case of the KLM, the first order term is just the 
mean- field result ^JrUa, while in the second order, one 
finds (Eq. (3.12) in ref. 39): 
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(30) 



< • • • >''^) means mean- field averaging, while the q- 
dependent spin operator is defined as usual, 



S^ = J2sre-^^^-''^ia = +,-,z) 



(31) 



SS^ is a short-hand notation: 



OOq — Dq \^q) 



(1) 



(32) 



In the following we are interested in the local self-energy 
^a{E) ^ j^J2k ^k<y{E) only, for which we find with (30) 
up to order J^ in the limit n ^ to be : 

-m,{m, + l)}GoiE) + O (J^) . 

(33) 

E. High-energy expansions 

For controlling unavoidable approximations, the spec- 
tral moments M^^ of the spectral density Sko- (E) 

Si^^E) = --ImGk^iE) (34) 

are of great importance: 

If dE-E^-S^AE). (35) 
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In principle, they can be calculated rigorously via the 
equivalent expression 
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There is a close connection between the spectral mo- 
ments and the high-energy behaviour of the Green func- 
tion: 
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(37) 



Because of the Dyson equation, 

EG^AE) ^h+ (e(k) + Sk.(£;)) G^AE) 
an analoguous expansion holds for the selfenergy: 
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The coefficients C^^- turn out to be simple functions of 
the moments up to order m+ 1: 
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(42) 



Using the definition (36), the moments of the KLM 
can be explicitly calculated by the use of the modcl- 
Hamiltonian (1). After tedious but straightforward ma- 
nipulations, one finds in the low- density limit (n ^ 0), 
for the first four moments: 



M <°j = 1 



K'J = f(k) - :^Jma 



(43) 



(44) 



K^J = e'(k) - Jm,e(k) + i J2(5(5 + i) _ ^^) (45) 



MltJ = e-^(k) - Ijm^e^ik) + ^J^A^ik) 



1 

+B{k) + 2e{k)iSiS + l)-m,)}+ (46) 
1 
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F. Interpolation formula 

We want to construct an approximate expression for 
the electronic selfenergy of the low- density KLM, which 
fulfills the zero-bandwidth limit (26) for all temperatures 
T and arbitrary coupling strengths J, as well as the exact 
T = 0-result (27) and (28) for arbitrary bandwidths and 
couplings. Furthermore, it should reproduce the correct 
high-energy (strong coupling) behaviour (39) and in ad- 
dition also the weak-coupling result (33). Guided by the 
non-trivial (n = Q,T — 0)-result (28), we start with the 
following ansatz for the local selfenergy: 



^AE) = -'-Jm^ + lj^ a.Go(i?-iJm.) 



4 1 - b^Go{E - Urn,) 



(52) 



Oct and feg. are at first unknown parameters. It is easy to 
recognize that this ansatz reproduces the exact limit (27) 
and (28) of ferromagnetic saturation, if T = 0, n = 0, 

flg- = (1 — Za-)S ; 5| = -J ; fe| arbitrary (53) 

and the zero-bandwidth limit (26), if 
e(k) ^ To Vk 



flcr = S{S + 1) — in„{mc + 1) 



(54) 



A(j(k) and B(k) are related to spin-correlation functions: 



«,j 
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hi 

Inserting these expressions into eqs. (40-42) we get for 
the first three selfenergy coefficients: 



^^J = -^rn^ 
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(49) 



Cl'J^-^J'iS{S + l)-m,{m, + l)) (50) 



CL'J = -^J'iM^) + B{k) e{k)ml) + 
+ Ij^{1 + m„){S{S + 1) - m^{m^ + 1)). 



(51) 



They determine the high-energy behaviour of the selfen- 
ergy (39). 



ba = b-a = -zJ 

We note that (54) agrees with (53) for T=0! 

By (52), we concentrate ourselves from the very be- 
ginning on the local part of the selfenergy. As already 
stated above, the wave-vector dependence of the selfen- 
ergy is mainly due to magnon energies 7iw(q) appearing 
at finite temperature in magnon emission and absorption 
processes by the band electron. However, the neglect of 
a direct Heisenberg exchange between the localized spins 
in the KLM can be interpreted as the h,uj{c\) -^ limit. 

We fix the parameters a^ and ba in the ansatz (52), 
by equating the high-energy expansion (39). For this 
purpose, we first develop (52) in terms of powers of the 
inverse energy. That requires the respective high-energy 
expression of the mean-field propagator Go{E — ^Jm^), 
which is exactly known: 
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it then follows: 
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The so-dcrivcd local selfenergy coefficients, 
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can be compared to the exact expressions following from 
(49)-(51): 



C(") = -\jm. 



(62) 



C« = ij2(5(5+l)-m,(m, + l)) (63) 



ci 



2) 



jj» I T„ 



w 



{S{S + 1) 
-m„{mcr + 1)) . 



(64) 



.(0) 



Cct is identically fulfilled. Agreement for the two other 

(65) 



coefficients is achieved by setting 

ftcr ~ S{S + 1) - ma{ma + 1) 



\^ 



b-n 



(66) 



These are the same expressions as found in (54) for the 
special zero-bandwidth limit. 

Inserting (65) and (66) into (52) yields a self-energy 
result, which is exact for T = (too- = ZaS) but arbi- 
trary bandwidths W and exchange couplings J. It fulfills 
the zero-bandwidth limit for all couplings J and all tem- 
peratures T. It obeys the high-energy behaviour which 
is important for the strong-coupling regime. Further- 
more, the comparison with (33) shows that the approach 
fits second-order perturbation theory, thus being reliable 
in the weak-coupling regime, too. We believe that (52) 
together with (65) and (66) represents a trustworthy ap- 
proach to the low-density self-energy of the Kondo-lattice 
model. In the next Section we present a numerical eval- 
uation. 



FIG. 1: Quasi-particle density of states as a function of energy 
for various values of magnetization. Full line for spin up and 
dotted line for spin down. J — 2, S — 7/2 and W = 1. 



III. RESULTS 

We have evaluated our theory for an sc lattice us- 
ing the respective Bloch-density of states (B-DOS) in 
tight-binding approximation [40]. The center of grav- 
ity To of the Bloch-band is chosen as energy zero. Fig.l 
shows the temperature-dependent quasiparticle density 
of states (Q-DOS) Pa{E) for a strongly coupled system 
(J = 2.0 eV, S = 7/2, 14^ = 1 eV^). The electronic 
spectrum gets its temperature-dependence exclusively 
through the local- moment magnetization m = \ma\ = 
(5*^)1, which must be considered as an external parame- 
ter. TO — 3.5 means T = QK (ferromagnetic saturation), 
while m = belongs to T = T^. The Q-DOS for each 
spin direction consists of two subbands separated by an 
energy of the order ^J{2S + 1). They originate from the 
two atomic levels Ei^ and E2a in the zero-bandwidth 
limit (16). 

A special case is the ferromagnetic saturation, for 
which the "f-spectrum consists only of the undeformed 
low-energy band {p^{E) = po{E + \JS)). The t 
—electron has no chance to exchange its spin with the 
perfectly aligned local-spin system. The spinflip terms 
in the exchange interaction (4) therefore do not work, 
only the Ising-like part (first term in (4)) takes care for a 
rigid shift of the excitation spectrum. The J, —spectrum 
is more complicated because a j-electron can, even at 
T — QK ^ exchange its spin with the ferromagnetically 
saturated spin system. One possibility is to emit a 
magnon, therewith reversing its own spin and becom- 
ing a I —electron. Such a spinfiip-excitation is of course 
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FIG. 2: Same as in figure 1 but with enlarged vertical scale 



possible only if there are | —states within reach on which 
the original \ —electron can land after the spinflip. That 
is the reason why the low-energy | — subband occupies 
the same energy region as the \ —band. 
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FIG. 3: Quasi-particle density of states(in the positive half of 
the frame) and imaginary part of the selfenergy(in the neg- 
ative half of the frame) as a function of energy for various 
values of magnetization. Full line for spin up and dotted line 
for spin down. J = 2, S' = 3/2 and W = \. 

The J, —electron has another possibility to exchange its 
spin with the ferromagnetically saturated moment sys- 



FIG. 4: Quasi-particle density of states as a function of energy 
for various values of magnetization. Full line for spin up and 
dotted line for spin down. J = 2, S = 3/2 and W = 1. 

For T > (to < 3.5) the "f —spectrum too, be- 
comes more structured because the localized spin sys- 
tem is no longer perfectly aligned. The system now con- 
tains magnons which can be absorbed by the ] —electron. 
Even polaron formation becomes possible. The spectral 
weight of the upper f — quasiparticle subband raises with 
increasing temperature, i.e., increasing magnon density. 
Fig.l illustrates that the temperature-dependence of the 
Q-DOS mainly affects the spectral weights of the sub- 
bands and not so much their positions. This is a typ- 
ical feature of the strong coupling regime JS » W . 
In such a situation, the band electron mobility is rather 
poor, it stays for a relatively long time at the same lat- 
tice site. The actual quantization axis is then the local- 
ized spin [S — 7/2), to which the electron can orient its 
spin parallel ("spin up" in the local frame) or antiparal- 
lel ("spin down" in the local frame). The excitation en- 
ergy for a parallel alignment roughly amounts to —-^JS, 
and for an antiferromagnetic alignment, to +2J{S + 1). 
The lower quasiparticle subband consists of states be- 
longing to the situation where the band electron appears 
in the local frame as "spin up" electron. This may hap- 
pen directly or after emitting/absorbing a magnon. In 
the upper subband the electron has entered the local 



frame as "spin down" electron. This is impossible for 
a 1" —electron at T = OK, when all localized spins are 
parallelly aligned {m = S). While the excitation ener- 
gies are almost temperature-independent, the probability 
for the electron to be in the local frame as a "spin up" 
or as a "spin down" particle strongly depends on tem- 
perature. That manifests itself in the spectral weight of 
the respective quasiparticle subband, which therefore is 
temperature- and spin-dependent. There remains a small 
probability that the band electron is not trapped by the 
localized spin, but rather propagates with high mobil- 
ity through the spin lattice. In such a case the effective 
quantization axis is no longer the local spin but rather the 
direction of the global magnetization (S^). Fig. 2 shows 
the Q-DOS for the same parameters as in Fig. 1 but 
on a finer scale. One recognizes two tiny satellites which 
emerge from the two main peaks with increasing temper- 
ature (decreasing magnetization m). The J,-satellite has 
a lower energy than the "f —satellite. This can be under- 
stood as follows: The original | —electron, for entering 
the low energy part of the spectrum, will predominantly 
do this by emitting a magnon thereby reversing its own 
spin. In case of being not trapped by a local spin , it 
then moves as a t —electron through the spin lattice. On 
the other hand, an original j —electron has to absorb 
a magnon in order to enter the high-energy part of the 
spectrum and propagating then as | —electron. With de- 
creasing magnetization the two satellites collapse mean 
field-like. In the strong coupling regime {JS >> W), 
pictured in Figs. 1 and 2, the satellites have only very 
small spectral weights, nevertheless representing interest- 
ing physics. 

The parameters used in Figs. 3, 4 and 5 (J = 
2 eV, W ~ 1 eV, S = 3/2) should be typical for the man- 
ganites. It is sometimes claimed [21,30] that because of 
the strong coupling JS, the itinerant electron (e^) spin is 
oriented at T = OK in any case parallel to the localized 
(t2g) spin. According to the exact m = S = 3/2— part 
of Fig. 3, this can strictly be ruled out for the FKLM. 
In the mentioned papers the assumption of full polar- 
ization is an artefact due to the restriction to classical 
spins (5* -^ oo)). The temperature-dependence of the Q- 
DOS is of course very similar to the S = 7/2 case in Fig. 
1. Even the satellites, which describe the "free" electron 
propagation after emitting/absorbing a magnon, do ap- 
pear (Fig. 4). However, because of the smaller distance 
between the two main peaks (« ^ J(2S' + 1)) the mean- 
field shift of the satellites is not so clearly visible as for 
the higher spin in Fig. 1. 

The imaginary part of the self-energy is directly re- 
lated to quasiparticle damping and lifetime, respectively. 
Fig. 3 demonstrates that the polaron states (upper part 
of the spectrum) represent quasiparticles with almost in- 
finite lifetimes since ImHaiE) is zero in this region . For 
T = OK, this is an exact result. In ferromagnetic satu- 
ration the whole t —spectrum consists of stable states. 
It turns out that, in the here discussed strong coupling 
regime, even for finite temperatures, only the states of 
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FIG. 5; Real and imaginary parts of the selfenergy as a func- 
tion of energy for different values of magnetization. For imag- 
inary part, full line for spin up and dotted line for spin down. 
For real part, dashed line for spin up and dash-dotted line for 
spin down. Imaginary part of the selfenergy is multiplied by 
a factor of 5 for better clarity. J = 2, S = 3/2 and W — 1. 



the mean-field satellites are getting finite lifetimes. The 
sharp peak of I'mT.^{E) falls always into the bandgap, 
which is provoked by a divergence of the real part of the 
selfenergy (Fig. 5). It has therefore no direct influence 
on the lifetime of quasiparticles. 

Up to now we have only discussed the FKLM in 
the strong coupling regime. As demonstrated in Sect. 
II. D, our interpolating approach is correct in the weak- 
coupling region, too. Fig. 6 shows, as an example, the 
Q-DOS for J = 0.2ey, W = leV and S = 3/2. The ten- 
dency to the two-subband structure can be recognized 
for weak couplings, too. The physical interpretation of 
the responsible elementary processes is the same as in 
the above-discussed strong coupling case. For T — 
all I —states represent stable quasiparticles, the respec- 
tive imaginary part of the selfenergy vanishes. With in- 
creasing demagnetization of the local moment system, 
IrnSij(E) becomes finite indicating finite lifetimes of 
t —quasiparticles due to magnon absorption, which is 
impossible at T = because of ferromagnetic saturation. 
Magnon emission by J, —electrons, however, is always 
possible. It should be pointed out that the upper part of 
the I —spectrum obviously consists, at low temperature, 
of stable polaron states. 

It is surprising that already very small couplings are 
sufhcient to create a pseudo gap in the quasiparticle 
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FIG. 6: Quasi-particle density of states(in the positive half of 
the frame) and imaginary part of the selfenergy(in the neg- 
ative half of the frame) as a function of energy for different 
values of magnetization. Full line for spin up and dotted line 
for spin down. Imaginary selfenergy is multiplied by a factor 
of 5 for better clarity. J = 0.2, S = 3/2 and W ^ 1. 
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FIG. 7: Quasi-particle density of states as a function of energy 
for different values of coupling constant J. Full line for spin 
up and dotted line for spin down, m = 1.5 (ferromagnetic 
saturation), S = 3/2 and W = 1. 



spectrum. According to Fig. 7, which shows the exact 
T = — Pa{E) for various exchange couplings J, the gap 
opens already for J ~ QAeV . Our rcsuhs for the weakly 
coupled FKLM are very similar to those presented in ref. 
10. 



IV. CONCLUSIONS 

We have presented an approach to the ferromagnetic 
Kondo-lattice model in the low-density limit {n -^ 0). 
The theory uses an interpolation formula for the elec- 
tronic selfenergy which fulfills a maximum number of lim- 
iting cases. It reproduces the non-trivial rigorous special 
case of a single electron in an otherwise empty conduc- 
tion band at T = (ferromagnetically saturated semicon- 
ductor) , and that for arbitrary bandwidths and coupling 
constants. It is exact in the zero-bandwidth limit for all 
temperatures and all exchange couplings. It obeys the 
high- energy expansion of the selfenergy, guaranteeing 
therewith the right strong-coupling behaviour, as well as 
perturbation theory of second order (oc J^) for the weak- 



tinct temperature-dependence in the electronic struc- 
ture, mainly due to spin exchange processes between 
the localized magnetic moments and itinerant band elec- 
trons. Magnon emission/absorption processes compete 
with polaron-like quasiparticle formation. These facts 
demonstrate that the assumption of classical spins (S — > 
oo), very often used for the simplified treatment of the 
model [30], suppresses just the essentials of Kondo-lattice 
model. The necessary extension of the presented theory 
has to include finite band occupations, which certainly 
requires additional approximations. The here developed 
n — > 0-approach can then serve as a weighty criterion for 
the correctness of the approach. 



V. ACKNOWLEDGEMENTS 

This work has been prepared as an India-Germany 
Partnership Project sponsored by the Volkswagen Foun- 
dation. 



1 W. Nohing, phys. stat. sol. (b) 96, 11 (1979) 
^ P. Wachter, in Handbook on the Physics and Chemistry of 
Rare Earths, edited by K. A. Gschneidner and L. Eyring ( 



Elsevier, Amsterdam, 1979), Vol. 2, p. 507 
^ G. Busch, P. Junod, P. Wachter , Phys. Rev. Lett. 12, 11 
(1964) 



11 



" J. Kossut, phys. stat. sol. (b) 78, 537 (1976) 

^ L. W. Roeland, G. J. Cock, F. A. MuUer, C. A. MoUman, 

K. A. M. Mc Ewen, R. C. Jordan, D. W. Jones, J. Phys. 

F 5 L 233 (1975) 
® S. Rex, V. Eyert, W. Nolting , J. Magn. Magn. Mat. 192, 

529 (1999) 
'^ Magnetism and Electronic correlations in Local Moment 

Systems: Rare Earth Elements and compounds edited by 

M. Donath, P. Dowben , W. Nolting (World Scientific, 

Singapore, 1998) 
® D. Li, J. Zhang, P. A. Dowben , M. Onellion, Phys. Rev. 

B 45, 7272, (1992) 
^ B. Kim, A. B. Andres, J. L. Erskine, K. J. Kim, B. N. 

Harmon, Phys. Rev. Lett. 68, 1931 (1992) 
^° W. Nohing, S. Rex, S. Mathi Jaya, J. Phys. Condens. Mat- 
ter 9, 1301 (1997) 
" S. Jin, T. H. Tiefel, M. Mc Cormack, R. A. Fastnacht, R. 

Ramesh, L. H. Chen, Science 264, 413 (1994) 
^^ A. P. Ramirez, J. Phys.: Condens. Matter 9, 8171 (1997) 
^^ W. Nolting, Quantentheorie des Magnetismus II, 

(Teubner-Verlag, Stuttgart, 1986), ch. 5.3.3. 
" C. Zener, Phys. Rev. 81, 440 (1951) 
^'^ S. Satpathy, Z. S. Popovic, F. R. Vukajlovic, Phys. Rev. 

Lett 76, 960 (1996) 
^^ W. E. Pickett, D. J. Singh, Phys. Rev. B 53, 1146 (1996) 
I'' D. J. Singh, W. E. Pickett, Phys. Rev. B 57, 88 (1998) 
i« J. H. Park, C. T. Chen, S. W. Cheong, W. Bao, G. Meigs, 

V. Chakaria, Y. U. Idzerda, Phys. Rev. Lett. 76, 4215 

(1996) 
^^ T. Saitoh, A. Sekiyama, K. Kobayashi, T. Mizokowa, A. 

Fujimori, D. D. Sarma, Y. Takeda, M. Takamo, Phys. 

Rev. B. 56, 8836 (1997) 



^° Y. Okimoto, T. Katsufuji, T. Ishikawa, A. Urushibara, T. 

Arima, Y. Tokura, Phys. Rev. Lett. 75, 109 (1995) 
^^ A. J. Millis, R. Miiller, B. I. Shraiman, Phys. Rev. B 54, 

5405 (1996) 
^^ E. Dagotto, S. Yunoki, A. L. Molvezzi, A. Morco, S. Cap- 

poni, D. Poilblanc, N. Furukawa, Phys. Rev. B 58, 6414 

(1998) 
2^ N. Furukawa, J. Phys. Soc. Jpn. 63, 3214 (1994) 
24 S. Doniach, PhysicaB 91, 231 (1977) 
'^^' H. von Lohneysen in ref. [7] 
2*5 P HiU, F. Willis, N. Ah , J. Phys.: Condens. Matter 4, 

5015 (1992) 
^■^ D. Gignoux, J. C. Gomez-Sal, Phys. Rev. B 30, 3967 

(1984) 
•^* N. Furukawa,in Physics of Manganites, edited by T. A. 

Kaplan and S. D. Mahanti, (Plenum, New York, 1999), p. 

1 
2^ A. J. Millis, P. B. Littlewood, B. I. Shraiman, Phys. Rev. 

Lett. 74, 5144 (1995) 
3° K. Held, D. VoUhardt, Phys. Rev. Lett. 84, 5168 (2000) 
^^ D. Meyer, C. Santos, W. Nolting, J. Phys.: Condens. 

Matter 13, 2531 (2001) 
^2 W. Nohing, M. Matlak, phys. stat.sol. (b) 123, 155 (1984) 
^3 B. S. Shastry, D. C. Mattis, Phys. Rev. B 24, 5340 (1981) 
'^^ S. R. Allan, D. M. Edwards, J. Phys. C 15, 2151 (1982) 
35 W. Nolting, U. Dubil, M. Matlak, J. Phys. C 18, 3687 

(1985) 
^^ H. Mori, Prog. Theor. Phys. 33, 432 (1965) 
^■^ H. Mori, Prog. Theor. Phys. 34, 399 (1966) 
3* G. Bulk, R. J. Jelitto, Phys. Lett. A 133, 231 (1988) 
3^ G. Bulk, R. J. Jelitto, Phys. Rev. B 41, 413 (1990) 
4° R. J. Jelitto, J. Phys. Chem. Solids 30, 609 (1969) 



